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Abstract. - A simple model for a superconductor with tunable critical current is studied 
theoretically. The model consists of a thin superconducting film with one vortex interacting 
with one magnetic dipole, whose magnetic moment is free to rotate, in the presence of a 
magnetic field applied parallel to the film surfaces. The pinning potential for the vortex is 
calculated exactly in the London limit. It is found that, due to the dipole freedom to rotate, 
the dependence of the pinning potential on the applied field is non-trivial, and allows both 
the spatial dependence and strength of the pinning potential to be changed by the field. As 
a consequence, the critical current can be tuned by the applied field. The critical current is 
obtained numerically as a function of the applied field. Order of magnitude changes in the 
critical current resulting from changes in the direction and magnitude of the applied field are 
reported, with discontinuous changes taking place in some cases. Possible application to vortices 
in low-T c superconducting films pinned by arrays of magnetic dots are briefly considered. 



The pinning of vortices in superconducting films by arrays of magnetic dipoles placed in 
the vicinity of the film is a topic that has received a great deal of attention lately. Most 
of the experimental [1] and theoretical [2-7] work carried out so far deals with arrays of 
permanent dipoles, that is, dipoles with magnetic moments fixed both in magnitude and 
direction. A related topic that has received little attention is vortex pinning by arrays of 
dipoles with magnetic moments free to rotate. The feasibility of fabricating such arrays has 
been demonstrated recently by Cowburn, et. al. [8]. These authors reported on the magnetic 
properties of arrays of nanomagnets made of Supermalloy, each nanomagnet being a thin 
circular disk of radius R. They found that for R ~ 50 — lOOnm the magnetic state of each 
nanomagnet is a single domain one with the magnetization parallel to the disk plane, and that 
the magnetization can be reoriented by small applied fields. One possible source of interest 
in vortex pinning by freely rotating dipoles is, as demonstrated in this paper, that the critical 
current may be tuned by an applied field. This paper studies theoretically in the London limit 
the interaction between one vortex in a thin superconducting film with one dipole, located 
outside the film, in the presence of a magnetic field parallel to the film surfaces. The magnetic 
dipole moment is assumed to be parallel to the film surfaces, to have constant magnitude and 
freedom to rotate. Tuning of the critical current is this model results because the interaction 
between the vortex and the dipole depends on the dipole orientation which, in turn, depends 
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on the applied field. Besides, in a thin film, a magnetic field parallel to the film surfaces has 
no effect on the vortex in the absence of the dipole. As shown here, this mechanism allows 
the pinning potential to be changed by the applied field over a wide range. When a transport 
current is applied to the film, the magnetic field created by it is parallel to the film surfaces 
and also contributes to the dipole orientation. This makes the pinning potential dependent 
on the transport current, and has important consequences for the critical current, as shown 
here. The main new results reported in this paper are: i) the exact analytic calculation of the 
pinning potential for one vortex interacting with a freely rotating dipole, and its dependence 
on the applied field and transport current, ii) the numerical calculation of the critical current 
for one vortex pinned by the dipole, and its dependence on the magnitude and direction of the 
applied field. This paper argues that these results are relevant for vortex pinning by arrays 
of nanomagnets, similar to those reported in Ref. [8], placed on top of superconducting films 
made of homogeneous materials, like most low-T c ones. The model is not applicable to layered 
high-T c superconducting films. 

The calculation of the pinning potential proceeds as follows. The superconductor film is 
assumed to be planar, with surfaces parallel to each other and to the x — y plane, isotropic, 
characterized by the penetration depth A, and of thickness d <C A. A vortex with vorticity 
q is located at position r, and the dipole is at ro = (0,0, zq > 0). The dipole moment m, 
has constant magnitude, m, is oriented parallel to the film surfaces, and is free to rotate in 
the x — y plane. An uniform magnetic field H is applied parallel to the film surfaces. The 
vortex-dipole system is shown in Fig. ^ The total energy in the London limit, neglecting 
pinning by random material defects, can be written as [7] 

E T = -m ■ (hi + H) + mH (1) 



where is the component parallel to the film surfaces of the field generated by the vortex at 
the dipole position. The energy Et does not include the vortex self-energy nor the interaction 
energy of the dipole with the field of the screening current generated by it in the film, because 
both are independent of the vortex position and dipole orientation. The constant mH is added 
for future convenience. The parallel component of the vortex field is given in the thin film 
limit (d < A) by [9] 
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This expression is exact for a thin film provided that r <C A = 2X 2 /d , which is the region of 
interest here. The total energy, Et, depends both on the vortex position r and on the dipole 




Fig. 1 - Superconducting film with one vortex at r, a magnetic dipole, m, at ro = (0,0, zo), and an 
applied magnetic field, H, parallel to the film surfaces. 

orientation. The pinning potential for the vortex at zero temperature, denoted by U vm , is 
the total energy for the equilibrium dipole orientation, that is, for m which minimizes Ex , 
with the vortex held fixed at r. Thus, according to Eq. QJ, the equilibrium m is parallel to 
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+ H, and the pinning potential is given by 

U vm = -m | + H | +mH . (3) 

Note that, by definition, U vm vanishes in the absence of a vortex. According to Eqs. © and 
(J5J, the spatial dependence of U vm is anisotropic. It depends both on r and on the angle 
between r and H. An important consequence of the dipole freedom to rotate is the non-trivial 
dependence of U vm on H obtained in Eq. ©. According to it H plays the role of a handle 
that controls the strength and spatial dependence of U vm , as will be discussed shortly. The 
scale for H in Eq. © is the vortex field, which is bound by b s ± < b s max = 0.3d/47rz x (0 O /A 2 ). 
It is convenient to use the following natural scales for physical quantities. Energy: eod, where 
e o = (</W47rA) 2 is the basic scale for energy/length of the superconductor. Magnetic moment: 
ftozo. Magnetic field: (f>o/\ 2 . 




Fig. 2 - Spatial dependence of the vortex pinning potential, U vm , (in units of eod) for q = 1, d = 
zo = 2£ , A = 10£ (£ =vortex core radius), m = 0.5<fioZo, and external field in the ^-direction, a) 
Permanent dipole. b) H = O.Ol^o/A 2 . c) H = O.Q2(j>o/\ 2 . d) H = 0. x and y in units of £. 

For H 3> b s rnax the dipole equilibrium orientation is parallel to H, and U vrn reduces to the 
pinning potential for a vortex interacting with a permanent dipole. Assuming that H is along 
the x-direction, U vm — ~mb xl which, according to Eq. J2J), coincides with the expression 
obtained in Refs. [3,4,6,7]. In this case U vm is anti-symmetric with respect to both an 
inversion of the vortex position (r — ► — r ==> U vrn — » —U vm ), and to a change in the sign 
of the vorticity ( q — * —q). For a vortex (q > 0), U vm has a minimum ( maximum ) on the 
x-axis at x = — (+)1.3zo, with minimum (maximum) value U vm = — (+)0.3 x ATreod(m/<poZo), 
as shown in Fig. 0a. In general, for H ^ the minimum of U vm occurs when b^_ is parallel 
to H, that is when the vortex (anti-vortex) is on the negative (positive) x-axis. In this case, 
according to Eq. U vm — —mb x . As a consequence, the minimum of U vm for H ^ is 
identical to that for a permanent dipole. However, the spatial dependence of U vm is strongly 
dependent on H, as shown in Fig.^Jfor some values of H < b s max (b s max — O.O240 O /A 2 for the 
parameters in Fig. EJl ■ For H = 0, U vm is given by U vm = —m | b^ |. In this case, according 
to Eq. J5J, U vm is the same for vortices and anti- vortices, has circular symmetry, and is 
attractive with a repulsive core, as shown in Fig. [3d. The minimum of U vm is degenerate 



4 



EUROPHYSICS LETTERS 



on a circle of radius r = 1.3zo , and has the same minimum value as a permanent dipole 
(U vrn = -0.3 x 47re o d(TO/0oZo))- 

Now the critical current, J c , for a single vortex with vorticity q = 1 is considered. The 
effect of a transport current density, J, applied to the film is twofold: it exerts on the vortex a 
force Fi = q(</>od/c)J x z and creates a field at the dipole position H; = (27rd/c)J x z, which 
adds to the external field and modifies the vortex pinning potential, because U vm is now given 
by Eq. with H replaced by the total field Hr = H + Hj. The critical current depends 
on the relative orientation of J and H. Here it is assumed that J is fixed in the positive 
y-direction, so that both Fl and Hj are along the positive ir-direction, and have magnitudes 
Hj = 2irdJ /c and Ft, = (fiodJ/c, and that H points in a direction that makes an angle a with 
the positive cc-axis, that is with F^. In this paper J c is obtained by solving numerically the 
equations of motion for the vortex. It is assumed that for J = the vortex is pinned at the 
absolute minimum of U vm , and that J increases very slowly with time. These assumptions 
ensure that the vortex follows the position of the minimum of U vm — Fl x as J increase, until 
J reaches a value for which the minimum becomes unstable, and the vortex depins. As J 
increases further, the vortex velocity also increases. The J c obtained here corresponds to J 
for which the vortex velocity reaches a small value chosen for numerical convenience. The 
obtained J c is slightly larger than the J for which the minimum becomes unstable. This is 
analogous to the voltage criterion in J c measurements. The values of J are, of course, limited 
to J < Jd, where Jd = c4>o / '(12\/37r 2 A 2 £) is the depairing current, £ being the vortex core 
radius. In the results reported next, regions where J c > Jd are discussed for the sake of 
completeness. Now there are two scales for H in U vm : b s maxl as discussed above, and Hj. The 
maximum Hj occurs for J — J c , and can be written as Hj c = 0.031d/£( J c /Jd)(4>o/^ 2 )- For 
d ~ zq ~ £, these two scales are comparable if J c ~ Jd- 




Fig. 3 - Single vortex (q = 1) critical current for A = 10£, d = zo = 2£: a) and b) J c vs. a; c) J c vs. 
H for constant a, indicated in the boxes; d) discontinuity field Hd vs m . Labels: m in units of (fioZo, 
H in units of 4>o/\ 2 . 

For H ^> (bfnax, Hj), U vm reduces to that for a permanent dipole oriented parallel to H, 
that is, with m making an angle a with the x-axis. In this case, U vm is independent of H 
and J and has a spatial dependence like that shown in Fig. 0a rotated by a with respect to 
the x-axis. For J = 0, the vortex is pinned at the absolute minimum of U vm , located at a 
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point in the x — y plane defined in polar coordinates, (p, 9), by (p = 1.3zq, 8 = a + tt) . The 
critical current depends on a and to, being a linear function of m, since U vm is linear in m. 
It is found that J c depends strongly on a, being largest for a = 0°, and decreasing smoothly 
with a, as shown in Fig. a) and b) ( curves labeled permanent dipole) . This results from 
the spatial dependence of U vm , as can be seen for a — 0°, 180°, where the critical current can 
be estimated analytically, because the vortex moves only along the x direction as J increases. 
The result is J c / Jd — im/<f>oZo for a = 0°, and J c / Jd — QAm/fozo for a = 180°. The origin 
of this tenfold difference can be seen in the plot of U vm shown Fig. 0a The driving force is 
parallel to the x-axis in Fig. 0a for a = 0°, and antiparallel for a = 180°. As can be seen in 
Fig. 0a, the slope of potential barrier is much steeper in the positive ^-direction than in the 
negative one. For other values of a the depinning process is more complicated, because the 
vortex motion as J increases is not confined to the direction of drive. 

For H comparable to or less than b s max and Hj, the equilibrium orientation of m is no 
longer fixed, and J c depends, besides on a and to, also on H. Typical results for A = 10. 0£ 
and d = zo = 2.0£ are shown in Fig. The J c vs. a curves are shown in Fig. 0a for 
to = 0.25^o^0) an d in Fig. 0b for to = O.50oz o , for characteristic values of H. In both cases 
the J c vs. a curves differ considerably from those for a permanent dipole for small H , being 
strongly dependent on H, and showing sharp changes in J c close to a — 180°, like those for 
to = 0.2500-zo, H = 0.001(/> /A 2 (Fig. 0a) and to = 0.5</> z , H = O.Ol^o/A 2 (Fig0b) . The 
curve labeled H = in Fig. 0a is the limit of the J c vs. a curve as H — > with a fixed. The 
strong dependence of J c on H is even more evident if J c is plotted as a function of H for fixed 
a, as shown in Fig. 0c for to = O.50o^o- In this case it is found that for a > 146.25° the J c vs. 
H curves have discontinuities at H = Hd, jumping from J c > Jd for H < Hd to J c ~ 0.2 Jd 
for H > Hd- For a < 146.25°, the dependence of J c on H is continuous, as illustrated by the 
curves for a = 135° and a = 90°. For a = 135°, J c undergoes a rapid change with H around 
Hd = 0.01A(f>o/ X 2 , whereas for a = 90° the change in J c with H is much slower. It is found 
that the J c vs. H curves have no discontinuities if to is smaller than a minimum value which 
depends on a. As shown in Fig. 0d, Hd vanishes at the minimum to, and increases above it 
essentially linearly with to. 




-1 1 2 3 0.0 0.5 1.0 1.5 

x/c J'J, 




Fig. 4 - Vortex positions with increasing J for m = O.50ozo, A = 10.0£, d — zo = 2.0£. A: Vortex 
initial position for J = 0. C: Positions where the vortex depins. a) and c) vortex trajectories. Dot 
indicates dipole location, b) and d) x and y coordinates vs. J corresponding to trajectories in a) and 
c). In d) top curves represent x/£, bottom curves y/£. Labels: H in units of 4>o/X 2 . 
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This complex behavior results from the dependence of U vm on J, through = H + Hj, 
as can be seen by examining how the position the minimum of U vm — Fi,x, which coincides 
with the vortex position, changes as J increases (Fig.Q). For a = 157.5°; H = 0.011</)o/A 2 < 
Hd, and a = 135°; H = 0.011</>o/A 2 , 0.013</>o/A 2 , when there are large enhancements in J c 
with respect to the permanent dipole value, the position of the minimum undergoes a large 
displacement, from the initial one on the right side of the dipole (A in Figs.0]a and0]c) to the 
final one, where the minimum becomes unstable, on the left side of the dipole (C in Figs. 0Ja 
and^c). This is accompanied by a flip in the direction of Hj- from near the negative x-axis at 
J = to one near the positive x-axis when the minimum becomes unstable. The enhancement 
in J c results because the vortex is effectively pinned by a permanent dipole oriented at a small 
angle with the positive x-axis. This can be seen for a — 157.5°; H = O.O110o/A 2 ( Fig. 0Jb), 
which shows that most of the vortex displacement from A to C takes place for < J < 0.5 J^. 
In this interval the direction of Ht rotates from 157.5° to 12° with the x-axis. When the vortex 
depins, at J = 1.35Jd, Ht points at 3° with the rr-axis and has magnitude Ht = 0.074</)o/A 2 . 
When there is little or no enhancement in J c (a — 157.5°; H = 0.0115</>o/A 2 , and a = 
135°; H — 0.015</>o/A 2 ) the position of the minimum undergoes only a small displacement, 
from A to B in Figs.0Ja and^Jc, and Ht points in a direction away from the x-axis. 

The reason for the discontinuous jumps in J c is related to way that the stability of the 
minimum of U vm — F^x changes as J increases. It is found that for H > Hd the minimum 
becomes unstable twice, whereas for H < Hd it becomes unstable only once. For H > Hd 
( a — 157.5°; H — 0.0115</>o/A 2 in Fig. 0Ja ) the minimum becomes unstable at B, where 
J = 0.25 3d- A stable minimum, not shown in Fig.0]a, appears again at a slightly larger value 
of J, and follows a trajectory close to the A-C curve. However, the vortex depins when the 
minimum becomes unstable for the first time at point B. For H = O.O1150o/A 2 > Hd in Fig. 
the minimum only becomes unstable once at point C. 

The J c results described above are believed to be representative of low-T c superconducting 
films. First, the particular set of parameters used, d *~ zq ~ £, are typical ones for supercon- 
ducting films with magnetic dots placed on top. For instance, in the experiments with arrays 
of magnetic dots with permanent magnetization placed on top of superconducting Nb films, 
reported in Ref. [10], d — 20nm ~ £. The magnetic dots are separated from the film by a thin 
protective layer of thickness ~ 20nm, so that the distance from the magnetic dipole to the 
film is zq ~ £. Second, since the dependence of J c /Jd on the model parameters d, zq, m, A, 
£, and H is, according to Eqs. ©, and only through the scaled variables d/zQ, m/4>oZo, 
and HX 2 /4>o, many superconducting film-dipole systems are equivalent. 

The London limit is valid for vortices in low-T c films. However, when a magnetic dipole is 
placed close to the film, it certainly breaks down if the dipole field destroys superconductivity 
locally in the film. Roughly speaking, London theory is valid as long as the maximum dipole 
field at the film is less than the upper critical field, that is, mj 'z^ < (f>o / (2tt^ 2 ) , or m/((f>QZo) < 
( z o I 'O 2 1 '■ For the values used in the above calculations (zo = 2£) this gives m/((f>oZo) < 0.64, 
which is larger than the values used in this paper. The London limit would be a better 
approximation if the present calculations were carried out for larger values of zo/C- However, 
the results for J c /Jd would be identical to those described above if m and d were scaled by the 
same factor as zo/£- For instance, if zq — > 2zq, J c /Jd would remain the same if d — > 2d and 
TTi — > 2m, but the upper limit of m/4>oZo for the validity of the London approximation would 
increase by a factor of 4. The present model also breaks down if m is sufficiently large to create 
vortices in the film. The threshold value of m for spontaneous vortex creation, estimated as 
777 ~ 0.7<fioZo using the results of Ref. [7], is larger tham m used here. 

The simple model discussed here is relevant to vortex pinning by arrays of magnetic dots, 
provided that: i) the dots are sufficiently far apart to neglect dipole-dipole interactions be- 
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tween them, ii) the number of vortices per dot is small enough, so that each dot pins at 
most one vortex, and the vortices are far enough apart to neglect vortex-vortex interactions. 
Unfortunately, there are no experimental results on vortex pinning by magnetic dots with 
freely rotating magnetic moments to compare the model predictions with. Instead, consider 
under which conditions the results described above apply to a system consisting of a typical 
array of nanomagnets reported in Ref. [8] on top of a thin superconducting film. Assuming 
that £ = 20nm, it follows that for d — z$ = 2£, A = 10£ (as above), d — zq — 40nm, 
A = 200nm, and <fio/\ 2 = 500G. The value m = 0.5</>oZo follows if the disk radius and 
thickness are chosen respectively as R ~ 50nm and t ~ lOnm, and the disk magnetization 
is taken as M ~ 10 2 /^s/(nm) 3 . If the distance between disks in the array is a ~ lfim, the 
dipole-dipole interaction energy , E^d ~ m 2 /a 3 , is small compared with the vortex pinning 
potential, U vm ~ —mb s max , since Edd/U vm ~ 10~ 2 . The values chosen for the disk radius 
and thickness, for the magnetization, and for the distance between disks are typical of those 
Ref. [8]. The results reported above (Fig. predict that for H < b s max = 12G, J c depends 
strongly on H , like in Fig. |31c, whereas for H > b s max — 12G, J c is that for a permanent 
dipole, and depends only on a. 

In conclusion then, this paper demonstrates that the critical current for a vortex in a 
thin superconducting film pinned by a freely rotating dipole can be tuned by a magnetic field 
applied parallel to the film surfaces. It is found that tuning takes place for a wide range of 
fields. For large fields, when the dipole moment is stuck in the field direction, the critical 
current changes continuously by one order of magnitude when the field is rotated 180°, from 
the direction parallel to the driving force to the direction opposite to it. For fields comparable 
to the vortex field the critical current is very sensitive to field variations, showing very rapid 
and even discontinuous changes by as much as one order of magnitude. It is suggested that 
the results apply to experiments on magnetic dot arrays on top of clean superconducting films. 

* * * 
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